##### vanilla simulation
source("~/evo-dispersal/People/Louis/HZDisp1D/HZDisp1DFunctions.R")

spX<-500

K<-20
lmda<-10
B<-1
strength<-1
n<-2*spX*K
ngens<-50
Dmean<-0
h2H<-1
h2D<-0.5
VPH<-0.05
VPD<-0.05

test<-init.inds(n, spX, Dmean, h2H, h2D, VPH, VPD)
test<-repro.disp(popmatrix=test, K=K, lambda=lmda, B=B, strength=strength, h2H=h2H,
	h2D=h2D, VPHI=VPHI, VPD=VPD, spX=spX)


out<-mother(n=n, spX=spX, K=K, lambda=lmda, B=B, strength=strength, ngens=ngens,
 Dmean=Dmean, h2H=h2H, h2D=h2D, VPH=VPH, VPD=VPD)


#X11()